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Abstract —  This  paper  focuses  on  the  average  path  length 
(APL)  of  BDD’s  for  switching  functions.  APL  is  a  metric  for  the 
time  it  takes  to  evaluate  the  function  by  a  computer  program. 
We  derive  the  APL  for  the  AND,  OR,  parity,  carry-out,  com¬ 
parison,  threshold  symmetric,  and  majority  functions.  We  also 
consider  the  average  of  the  APL  for  various  classes  of  functions, 
including  symmetric,  threshold  symmetric,  and  unate  cascade. 
For  symmetric  functions,  we  show  the  average  APL  is  close 
to  the  maximum  path  length,  n,  the  number  of  variables.  We 
show  there  are  exactly  two  functions,  the  parity  functions,  that 
achieve  the  upper  bound,  n,  on  the  APL  for  BDD’s  over  all 
functions  dependent  on  n  variables.  All  other  functions  have 
an  APL  strictly  less  than  n.  We  show  that  the  APL  of  BDD’s 
over  all  functions  dependent  on  n  variables  is  bounded  below 
by  2  —  2T11_1  ■  The  set  of  functions  that  achieves  this  small  value 
is  uniquely  the  set  of  unate  cascade  realizable  functions.  We  also 
show  that  the  APL  for  benchmark  functions  is  typically  much 
less  than  for  random  functions. 


I.  Introduction 

Within  the  past  25  years,  considerable  research  has  been 
devoted  to  minimizing  the  number  of  nodes  in  binary  deci¬ 
sion  diagrams  (BDD’s).  This  problem  is  motivated  by  the 
fact  that  the  amount  of  memory  needed  to  store  a  function  as 
a  BDD  is  directly  proportional  to  the  number  of  nodes.  How¬ 
ever,  there  is  another  cost,  the  cost  of  evaluating  the  function 
for  some  combination  of  variable  values.  This  is  the  subject 
of  this  paper. 

We  assume  that,  for  a  given  function  /,  all  2"  assignments 
of  values  to  the  n  variables  are  equally  likely.  Thus,  the  aver¬ 
age  path  length  (APL)  in  the  BDD  of  /  can  be  computed  by 
summing  the  path  lengths  over  all  2n  assignments  of  values 
to  the  n  variables  and  dividing  by  2n.  We  are  also  interested 
in  the  distribution  of  path  lengths,  namely  the  number  paths 
for  each  path  length  value.  Therefore,  we  can  determine  the 
extent  to  which  pathological  cases  occur.  To  illustrate,  con¬ 
sider  the  BDD  of  the  carry-out  function  of  a  6-bit  ripple  carry 
adder,  where  n  =  26  +  1.  That  is,  we  assume  each  bit  to  be 
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(a)  Carry-out  Cb  of  6-bit  ripple  carry  adder  circuit. 


(b)  BDD  for  Cb  with  MSB  at  top.  (c)  BDD  for  Cb  with  LSB  at  top. 
Fig.  1.  Two  BDD's  for  the  carry-out  function  of  a  6-bit  ripple  carry  adder. 


added  contributes  two  input  variables,  and  there  is  one  carry- 
in  variable.  Fig.  la  shows  the  ripple-carry  adder  circuit  that 
produces  this  function.  Fig.  lb  shows  the  BDD  in  which 
the  variables  are  in  descending  order  (most  significant  bits  at 
the  top),  and  Fig.  lc  shows  the  BDD  in  which  the  variables 
are  in  ascending  order  (least  significant  bits  at  the  top).  Note 
that  a  dotted  edge  corresponds  to  Xi  =  0  and  a  solid  line  cor¬ 
responds  to  Xi  =  1.  Although  both  orderings  yield  BDD’s 
with  the  same  number  of  nodes,  the  first  has  more  shorter 
paths,  as  shown  by  the  distributions  in  Fig.  2  for  n  =  33. 
The  sawtooth  pattern  associated  with  the  most  significant  bit 
(MSB)  at  the  top  occurs  because  there  are  no  paths  of  odd 
length.  In  this  distribution,  most  paths  have  lengths  less  than 
10.  The  distribution  associated  with  the  least  significant  bit 
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Fig.  2.  Distributions  of  path  lengths  for  the  carry-out  function. 

(LSB)  at  the  top  is  binomial  and  centered  around  path  length 
25.  Thus,  the  APL  for  the  LSB  at  the  top  is  25,  which  is  sig¬ 
nificantly  more  than  the  APL  for  the  MSB  at  the  top,  4.0.  In 
pass  transistor  implementations  [1]  [2]  of  BDDs,  significant 
differences  in  APL  result  in  significant  differences  in  delay 
in  logic  circuits. 

The  path  length  in  decision  diagrams  is  important  in 
databases  [3],  pattern  recognition  [4],  and  analysis  of  algo¬ 
rithms;  for  example,  decision  trees  have  been  used  to  deter¬ 
mine  lower  bounds  on  the  complexity  of  sorting  algorithms 
[5].  [6]  is  an  easily  understood  (but  old)  survey.  Recently, 
Liu,  Wang,  Hwang,  and  Liu  [7]  show  three  heuristics  for 
finding  the  minimal  APL  in  BDD’s.  They  show  that  the  or¬ 
dering  of  variables  that  achieves  1 .  the  minimum  path  length 
and  2.  the  minimum  size  (number  of  nodes)  is  different  for 
all  of  the  20  benchmark  functions  they  considered. 

Our  results  extend  the  results  of  Brandman,  Orlitsky,  and 
Hennessey  [8]  in  which  lower  bounds  on  the  APL  are  com¬ 
puted.  In  our  paper,  we  compute  exact  values  for  the  APL, 
as  well  as  distributions  of  path  lengths.  We  show  that  some 
distributions,  like  that  of  the  AND  function,  extend  over  the 
whole  range  of  possible  path  lengths,  while  others,  like  that 
of  the  parity  functions,  are  very  narrow. 

Because  of  space  limitations,  we  do  not  include  proofs. 
For  a  complete  version  of  this  paper,  see  [9] 

II.  Individual  Functions 

In  this  section,  we  consider  the  path  lengths  in  BDD’s  as¬ 
sociated  with  individual  functions.  Our  approach  is  to  use 
generating  functions  to  track  the  number  of  paths  of  various 
lengths. 

A.  AND  and  OR  Function 

Fig.  3a  shows  the  BDD  of  an  rc-variable  AND  function. 
The  generating  function  for  the  distribution  of  path  lengths  in 
the  (trivial)  BDD  of  the  AND  function  on  one  variable  is  2 z1. 
That  is,  there  are  two  assignments  of  variables  (x  \  =  0  and 
X\  =  1)  each  corresponding  to  a  path  of  length  1.  The  AND 
function  of  two  variables  has  a  path  length  distribution  de¬ 
scribed  by  2 z1  +  2 z2.  In  this  case,  two  assignments  of  values 


Fig.  3.  BDD’s  for  the  AND  and  parity  functions. 

cause  the  AND  to  be  evaluated  as  0  through  a  path  of  length 
1  ( X1X2  =  00  and  01)  and  two  assignments  (X1X2  =  10  and 
11)  cause  the  function  to  be  evaluated  as  0  or  1  through  a 
path  of  length  2.  For  the  n-variable  AND,  the  path  length 
distribution,  PLD(AND(n))  is 


PLD(AND(n ))  =  2n~1z  +  2n~2z2  +  ■  ■  ■  +  22zn~2  + 
21zn~1  + 2  •  2° zn 


Dividing  this  by  2"  yields,  FPLD{AND{n)),  a  generat- 
ing  function  for  the  fraction  of  paths  of  length  l,  2,  etc.. 

1  —  zn+1  zn 

FPLD(AND(n))  =  ^  -  L  0) 

The  APL  for  the  AND  function,  APL{AND{n))  is 
calculated  by  forming  the  weighted  sum  of  path  lengths 
and  dividing  by  2n,  the  number  of  assignments  of  values 
to  the  n-variables.  This  can  be  done  by  differentiating 
PLD(AND(n))  with  respect  to  z,  multiplying  by  z,  and 
setting  z  =  1.  For  example,  from  the  discussion  above, 
PLD(AND(2))  =  2 z1  +  2 z2.  Differentiating  with  respect 
to  2  and  multiplying  by  z  yields  2Z1  +  4 z2  which  has  the 
effect  of  generating  the  weighted  sum  associated  with  path 
of  lengths  1  (2Z1)  and  2  (4z2).  Setting  z  =  1  forms  the  sum 
over  all  weighted  values,  which  is  6  in  this  case.  Dividing 
by  22,  the  number  of  assignments  of  values  to  two  variables, 
yields  APL{AN D{ 2))  =  1.5.  Performing  these  steps  on  (1) 
yields  a  result  by  Breithart  and  Gal  [10] 

Lemma  2.1. 

APL{AND(n ))  =  APL(OR{n))  =  2  -  — 

The  above  lemma  also  applies  to  the  OR  function,  because 
its  BDD  is  isomorphic  to  that  of  the  AND  function. 


B.  Parity  Function 


Fig.  3b  shows  the  BDD  for  a  parity  function  (the  exclusive 
OR).  Here,  all  paths  from  the  root  node  to  a  terminal  node 
include  a  node  associated  with  each  variable.  Since  there  are 
no  ’’short-cuts”,  the  length  of  all  paths  is  n,  and  the  distribu¬ 
tion  of  nodes  is  2 nzn. 

Lemma  2.2. 

APL(PARITY(n ))  =  n. 

C.  Functions  With  the  Largest  and  Smallest  APL 

In  a  function  whose  APL  is  exactly  n,  all  paths  must  have 
maximum  length.  Functions  with  this  property  are  rare,  as 
shown  below. 


Lemma  2.3.  The  parity  functions  are  uniquely  those  func¬ 
tions  whose  BDD’s  have  the  largest  APL  (n)  among  all  n- 
variable  functions. 

The  AND  and  the  OR  functions  also  have  a  special  dis¬ 
tinction. 

Definition  2.1.  f  is  a  unate  cascade  realizable  function  if  f 
can  be  represented  as 

f{x l,x2,  ...Xn)  =  •  ■  «-l0n-lO)  •  •  •)> 

where  x*  is  either  x-i  or  xt  and  <C>i  is  either  the  OR  V  or  AND 
A  function. 

For  example,  if  all  variables  occur  uncomplemented  and 
all  operations  <),:  are  AND  (respectively,  OR),  the  resulting 
function  is  the  AND  (respectively,  OR)  of  all  variables.  The 
significance  of  a  unate  cascade  realizable  function  is  that  its 
BDD  has  the  property  that  all  nodes  have  at  least  one  edge 
from  that  node  to  a  terminal  node.  Further,  for  any  such 
BDD,  we  can  find  a  unate  cascade  realizable  function  cor¬ 
responding  to  it. 

Lemma  2.4.  The  unate  cascade  realizable  functions  are 
uniquely  those  functions  whose  BDD’s  have  the  smallest 
(2  —  2„1_1 )  APL  among  all  n-variable  functions. 

D.  Carry-out  and  Comparison  Functions 

In  this  section,  we  consider  the  carry-out,  X  =  Y  (equal 
to),  and  X  >  Y  (greater  than  or  equal  to)  functions.  Fig. 
4  shows  the  BDD’s  for  the  X  >  Y  and  X  =  Y  function, 
where  the  MSB  is  at  the  top.  Fig.  1  shows  the  BDD  for  the 
carry-out  function. 

In  the  case  of  the  X  =  Y  function,  the  position  of  the  pair 
( XiUi )  in  the  ordering  is  irrelevant.  However,  for  the  carry¬ 
out  and  X  >  Y  functions,  the  smallest  APL  occurs  when  the 
MSB  is  at  the  top.  For  example,  in  the  X  >  Y  function,  the 
MSBs  are  always  needed  to  determine  whether  the  output  is 
1  or  0.  Suppose  that,  the  lowest  MSB  from  the  top  variable 


Fig.  4.  BDD’s  for  the  X  =  Y  and  X  >  Y  functions. 

(either  Xf,  or  yf)  is  the  fc-th  variable  from  the  top,  and  assume 
that  k  >  2.  The  two  MSBs  are  needed  to  determine  the 
function  value.  Further,  the  other  (k  —  2  >  0)  variables  are 
needed  if  Xb  =  iji,.  It  follows  that  all  paths  from  the  root  node 
to  a  constant  node  have  length  k.  We  can  state. 

Lemma  2.5.  A  BDD  for  the  carry-out  and  X  >  Y  functions 
of  smallest  APL  is  achieved  by  placing  the  MSBs  at  the  top, 
the  next  most  significant  bits  next,  etc.. 

As  in  the  previous  examples,  we  derive  a  generating  func¬ 
tion  for  the  distribution  of  paths  to  path  lengths,  and,  from 
this,  derive  the  APL.  The  results  are  shown  in  Table  I.  It  is 
interesting  that  all  three  functions  have  an  APL  that  is  a  con¬ 
stant  4  for  large  n.  Except  for  the  bottom,  the  BDD’s  struc¬ 
ture  is  the  same  for  all  functions.  Indeed,  the  difference  in  the 
values  is  explained  entirely  by  the  structure  at  the  bottom. 

TABLE  I 

Distribution  of  path  length  and  APL  for  comparison 

FUNCTIONS 


Function 

Distribution 

APL 

Carry-out 

X  =  Y 

X  >  Y 

22bz2(l  -  227)/(l  -  ^)  +  2b+1z2b+1 
(22b~1z2  -  2bz2b)/(l  ~  4)  +  2 b+1z2b 
(2 2b~lz2  -  2b~1z2b) / (1  -  )  -  2 bz2b~l 

4^  4^ 

1  1 

E.  Symmetric  Threshold  Function 

A  symmetric  function  has  the  property  that  permuting  any 
of  the  variables  leaves  the  function  unchanged.  A  symmetric 
threshold  function  is  a  symmetric  function  that  is  1  iff  t,  of  the 
input  variables  are  1,  for  0  <  t  <  n.  The  majority  function 


Fig.  5.  BDD  of  the  5-variable  majority  function 


Lemma  2.10. 


is  a  special  case  of  a  symmetric  threshold  function,  where 
n  =  2t  —  1  and  t  =  1,  2,  ...  . 

For  this  class  of  functions,  we  can  state 

Lemma  2.6.  The  generating  function  for  the  distribution  of 
the  fraction  of  paths  with  various  path  lengths  for  the  sym¬ 
metric  threshold  function  with  threshold  t  is 


FPLD(S.THRES{n,t))  = 


n 


E 


j~  1 
k-  1 


j~  1 

n  —  k 


where  k  =  min{£,  n  —  t  +  1}. 

Lemma  2.7. 


fc  ("-A 

APL(STTHRES(n,t))  =  2fc  -  £  j. 

i=i 


where  k  =  min{f,  n  —  f  +  1}, 
and  when  n  is  large,  we  can  write 

Lemma  2.8. 


APL(S-THRES(n,  t))  ~  2  min{f,  (n  —  f  +  1)}, 
where  A(n)  ~  means  lim,^^  =  1. 


F  Majority  Function 

A  majority  function  is  a  special  case  of  a  symmetric 
threshold  function.  Specifically,  n  =  2t  —  1,  where  t  is  the 
threshold.  Fig.  5  shows  the  BDD  of  the  majority  function 
on  5  variables.  In  the  case  of  a  symmetric  threshold  func¬ 
tion,  we  assumed  that  the  threshold  t  stayed  fixed  as  n,  the 
number  of  variables,  increased.  In  the  case  of  the  majority 
function,  we  assume  that  t  increases  as  n  increases  (so  that 
the  relation  n  =  2f  —  1  always  holds).  We  have 

Lemma  2.9. 

FPLD(MAJ(n))  =  2 V  £  (|)’ .  (2) 

where  n  =  2t  —  1,  for  t  =  1,2,.... 


APL(MAJ(n ))  =  n  -  +  1,  (3) 

where  n  is  odd. 

The  second  term  in  (3)  contains  the  factor  (fji),  which  is 
„_i7,,;,+1 , .  The  factorials  can  be  approximated  using  Ster- 

~/T~ T~  • 

ling’s  approximation,  yielding 

Lemma  2.11. 


APL(MAJ(n)) 


n  — 


+  1. 


Here,  we  retain  the  term  +1  because  it  is  significant  for  small 
values  of  n.  Indeed,  this  approximation  is  quite  accurate  for 
small  values  of  n.  For  example,  for  n  =  3,  the  exact  expres¬ 
sion  and  the  approximation  differ  by  less  than  5%,  while  for 
n  =  15,  the  difference  is  only  0.4%. 


G.  Comparison  of  Distributions 

Fig.  6  shows  the  distributions  of  path  lengths  for  the  AND, 
majority,  and  parity  functions  on  n  =  33  variables.  The  path 
lengths  in  the  BDD  of  the  AND  function  range  from  1  to  33, 
while  the  path  length  distribution  of  the  parity  function  is  a 
single  point,  at  n  =  33.  In  between,  is  the  majority  path 
length  distribution. 


III.  APL  for  Sets  of  Functions 

In  this  section,  we  consider  the  APL  over  sets  of  functions. 
We  assume  that  each  function  contributes  equally  to  this  ’’av¬ 
erage  of  averages”.  We  consider  four  sets  of  functions,  all 
functions,  all  (totally)  symmetric  functions,  all  symmetric 
threshold  functions,  and  all  unate  cascade  realizable  func¬ 
tions. 

Definition  3.1.  APLs(n)  denotes  the  average  path  length 
of  BDD ’s  for  n-variable  functions  in  the  set  S.  That  is, 
APLs(n )  =  APL(fi)  ,  where  S  denotes  a  set 

of  functions. 


TABLE  II 

APLau(n )  for  3  <  n  <  16  obtained  by  averaging  10,000 

SAMPLES. 


3  f2. 187500 

4 1  f3. 183594 1 

FT[j  4.226113 1 

6  5.203656 

7  6.195206 

8  7.188487 

9  8.186332 

10  9.184930 

11  10.184335 

12  11.183719 

13  12.183499 

14  13.183635 

15  14.183664 

1 16  |  15.183647  | 

f  Exact  value  obtained  by  enumerating  all  functions. 

A.  Set  of  All  Functions 

In  this  section,  we  determine  an  upper  bound  on  the  APL 
for  all  functions  on  n  variables.  Let  APLau(ji,  II)  be  the 
APL  over  all  n-variable  functions  for  fixed  ordering  II  of  the 
variables. 

Lemma  3.1. 

APLau(ti,  II)  =  APLAii(n  —  1,11)  +  1  —  22n~1 ' 

Theorem  3.1. 

APLAU(n)  <(n-  1)  +  0.183578. 

Table  II  shows  the  APL  for  randomly  generated  functions 
on  n-variables,  for  3  <  n  <  17.  For  n  =  3  and  n  =  4,  exact 
values  are  obtained  by  enumerating  all  functions.  For  each 
value  of  n  in  the  range  5  <  n  <  14,  10,000  samples  were 
generated. 

Note  that  this  data  was  generated  by  choosing  some  fixed 
ordering  of  the  variables,  and  then  computing  the  APL 
through  the  sample  function.  This  data  matches  closely  the 
upper  bound. 

B.  Set  of  All  Functions  With  a  Specified  Number  of 
Minterms 

We  now  consider  how  the  number  of  true  minterms  in  a 
function  affects  its  APL.  That  is,  instead  of  considering  all 
n-variable  functions  as  a  single  set,  we  divide  this  set  into 
subsets,  according  to  the  number  true  minterms.  Then,  we 
seek  the  APL  of  all  functions  in  each  subset.  Table  III  shows 
APLau(n ,  II)  the  average  APL  of  BDD’s  for  n-variable 


functions  with  k  true  minterms.  Each  data  value  is  the  av¬ 
erage  of  10,000  random  samples.  The  APL  values  are  sym¬ 
metric  about  the  middle  k  value,  since  the  BDD  of  a  function 
with  k  true  minterms  is  isomorphic  to  the  BDD  of  its  com¬ 
plement,  which  has  n  —  k  true  minterms.  We  have  omitted 
this  mirror-image  data. 

table  in 

APLaii(n,  k ,  LI)  FOR  n  =  6,  8,  and  10  VARIABLES. 


k/  2" 

APLaU(n,  fc,II) 
n  =  6  n  =  8  n=10 

0.0000 

0.000000 

0.000000 

0.000000 

0.0625 

3.364556 

5.246588 

7.218377 

0.1250 

4.152775 

6.097209 

8.082283 

0.1875 

4.578462 

6.540764 

8.533054 

0.2500 

4.844106 

6.816895 

8.810673 

0.3125 

5.013888 

6.995677 

8.989050 

0.3750 

5.121319 

7.107553 

9.101985 

0.4375 

5.184331 

7.169411 

9.164830 

0.5000 

5.203656 

7.188487 

9.184930 

C.  Set  of  All  Symmetric  Functions 

Moret,  Thomason,  and  Gonzalez  [11]  show  that  a  BDD 
of  any  symmetric  functions  dependent  on  n  variables  has  at 
least  one  path  that  is  the  longest  possible,  n.  In  this  section, 
we  extend  this  by  deriving  the  average  APL  of  BDD’s  of 
symmetric  functions. 

Lemma  3.2. 

AP  L  gym  (p)  =  APLgym(n  —  1)  +  1  —  — 

Theorem  3.2. 

APLSym{n)  =  n  -  1  + 

D.  Set  of  All  Symmetric  Threshold  Functions 

A  symmetric  threshold  function  is  a  symmetric  function 
that  is  1  iff  k  or  more  of  the  variables  are  1.  For  example,  the 
OR  and  AND  functions  on  n  variables  are  symmetric  thresh¬ 
old  functions  with  k  =  1  and  k  =  n,  respectively.  We  have 

Theorem  3.3. 

o 

AP L gymThres  (it)  =  y  2)  ' 

IV.  Benchmark  Functions 

The  APL  of  a  randomly  generated  function  of  n  variables 
is  near  n.  It  is  interesting  to  consider  the  APL  for  benchmark 


Fig.  7.  Graph  of  APL  versus  n  for  4352  benchmark  functions. 

functions.  To  investigate  this,  we  considered  189  benchmark 
functions  in  [12],  For  each  output  f3  of  a  multiple -output 
function,  we  minimized  the  number  of  nodes  in  the  BDD  for 
fj  (using  an  exhaustive  method  for  n  <  16  and  a  heuristic 
method  for  16  <  n)  and  computed  the  resulting  APL.  In  all, 
we  minimized  4352  single-output  functions.  For  each  n,  we 
collected  the  functions  that  depend  on  n  variables,  and  ob¬ 
tained  the  average  of  their  APL.  To  compare  among  various 
values  of  n,  we  divided  this  average  by  n. 

The  result  is  shown  in  the  graph  of  Fig.  7,  in  which  the 
value  of  is  plotted  verses  n.  Each  x  in  this  figure  corre¬ 
sponds  to  the  average  of  all  n-variable  benchmark  functions 
normalized  by  dividing  by  n.  Note  that  APL  is  much  smaller 
than  1.0  except  for  small  values  of  n  and  n  =  41.  The  out¬ 
lier  at  n  =  41  is  due  to  function  C499,  used  for  single  error 
correction,  which  has  32  outputs  on  41  variables,  each  with 
an  APL  of  approximately  25.  There  is  a  tendency  for  ' —  -  to 
decrease  as  n  increases.  This  is  to  be  compared  to  randomly 
generated  functions,  where  APL  tends  to  increase  as  n  in¬ 
creases.  This  shows  that  the  benchmark  functions  have  quite 
different  properties  than  randomly  generated  ones.  Indeed, 
since  we  minimized  the  node  count  and  not  the  APL,  we 
would  expect  an  even  greater  difference  had  we  minimized 
the  APL. 


V.  Concluding  Remarks 

In  this  paper,  we  consider  the  average  path  length  or  APL 
of  BDD’s,  as  measured  by  the  average  number  of  edges  tra¬ 
versed  from  the  root  node  to  a  terminal  node.  We  have  shown 
that  the  APL,  for  some  functions,  such  as  the  carry-out,  is 
strongly  dependent  on  the  ordering  of  the  variables.  We  have 
derived  the  APL  for  various  classes  of  functions.  This  is  sum¬ 
marized  in  Table  IV 


TABLE  IV 

Average  path  length  in  BDD's  of  sets  of  functions  for  large 

n. 


Function 

APL 

PARITY 

n* 

AND  and  OR 

o  1  * 

z  2"-1 

Symmetric  Threshold  k 

2k 

MAJ  (Majority) 

n-  \[m 

V  7T 

BTREE  (Balanced  Tree) 

log2 

Carry-out 

4-^ 

All 

<  n  —  0.816 

All  Symmetric 

n  -  1  +  ± 

All  Symmetric  Threshold 

n  +n  * 
2(ra+2) 

All  Unate  Cascade 

O  1  :jc 

Z  2^  —  1 

*  Exact  value. 
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